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Abstract 


/^-synthesis is a widely used controller design paradigm that utilizes the structure in 
the uncertainty associated with the plant. However the controller thus produced is 
typically of very high order. This dissertation is concerned with the approximation 
of the controller such that there is no degradation in the closed loop performance. A 
new proof has been proposed for the additive perturbation reduction technique with 
sufficient conditions to guarantee the closed loop structured singular value to remain 
less than unity. This proof provides a rigorous basis for controller reduction in the 
/i framework. A new coprime reduction technique has been proposed that makes 
unstable controller reduction possible with sufficient conditions to guarantee preser- 
vation of closed loop performance. The coprime factor perturbations to the controller 
have been shown to have a block diagonal structure. The proposed algorithms have 
been tested on a widely studied benchmark HIMAT aircraft and have been found 
to work satisfactorily producing more than 50% reduction in the controller order 
without optimization. Lastly, a full ^ synthesis design has been performed on the 
experimental data available from a flexible launch vehicle. The details of the design 
with the selected performance weights and the simulation results are presented. The 
reduction algorithms are also tested on this practical example and found to produce 
considerable reduction in the controller order while preserving performance. 
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Introduction 


The first step in controlling a system is to identify the system. Through identifi- 
cation a set of mathematical equations is derived that should accurately describe 
the dynamics of the system. However physical systems are very high dimensional, 
and an inherent tradeoff exists between the accuracy of the representation and the 
complexity of the equations chosen to represent it. In this thesis, we consider only 
systems that can be represented by finite dimensional linear time invariant differen- 
tial equations. Thus the tradeoff mentioned above reduces to a minimization in the 
order of the differential equations that have been used to model the system without 
sacrificing the accuracy of such a representation. This so called model approxima- 
tion problem for FDLTI systems has been extensively studied. One of the widely 
practised methods in this regard is to define a set of low dimensional systems, which, 
as a whole encompasses the complex behaviour of the real system. Thus this entire 
set is considered as a representation of the the original “uncertain” system. In other 
words this modelling procedure decreases complexity by introducing uncertainty into 
the system. The most common procedure to define this set of plants is to define a 
nominal plant and then consider a set of perturbations about this nominal plant to 
define the entire set. Depending on the structure of this set of perturbations different 
methods in describing the uncertainty have been evolved. This concept of uncertainty 
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is fui'ther extended so as to encompass other unknown or difiicult-to-model parame- 
ters like nonlinearity and time variation etc. Under such a representation the control 
objective becomes designing a controller that provides stability and performance to 
all plants belonging to that set. Thus it is claimed that such a controller will be able 
to provide stability and performance when connected to the physical plant. Two of 
the most widely studied methods of controller design of such uncertain systems are 
the Tioo-synthesis and the /.t-synthesis procedures. (See Francis [7], Doyle et al [5] 
and Packard and Doyle [24] for introductions to 77oo and fi theory.) 

Unfortunately the controllers designed by this paradigms are usually of very high 
order (typically equal to or greater than the order of the nominal plant). Practical 
controllers must be simple, linear and of low order. High complexity controllers are 
not only difficult to understand but they may turn out impossible to implement 
in hardware and software. Also from the point of integrity and reliability, simple 
controllers are preferable. In this dissertation we are chiefly concerned with the 
reduction of the state dimension of controller designed in a framework such that the 
closed loop structured singular value, with respect to the given uncertainty structure, 
remains less than unity. This guarantees that there shall be no degradation in the 
Tioo performance of the closed loop system. 

The design of reduced order controllers for high order plants may be divided into 
three categories [1] : direct design via optimization of controller parameters; model 
reduction followed by the design of a low order controller using the reduced order 
model; and design of a high order controller from the high order model, followed by 
reduction of the high order controller. The first method is evidently most desirable 
and much work has been done in this regard [14], [15], [16], [21]. However the problem 
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remains mathematically unsolved. The second method is intuitively appealing but it 
does not provide any a priori guarantee that the controller that is working for the low 
order model shall work for the actual system. In this dissertation we concentrate on 
the third method. 

The main contribution of this thesis are as follows: a new proof for the Kavranoglu’s 
additive reduction technique [17] has been provided and the weights for two sided re- 
duction have been derived. A new coprime factor reduction scheme has been proposed 
which guarantees closed loop stabihty and performance with structured imcertainty. 
The coprime factor perturbations to the controller have been found to have a block 
diagonal structure thus improving the reduction algorithm. In Chapter 4 a complete 
yu-synthesis design is performed for the experimental data available from a flexible 
launch vehicle. 

The organization of the thesis is as follows: 

Chapter 1 In this chapter a brief description of the mathematical tools needed in 
the analysis of the Controller Reduction problem are discussed. Concise outflnes 
of the generic state space system, the Hardy Spaces, Internal Stability, Linear 
Fractional Transforms, Coprime and Spectral Factorizations, the Standard Tioo 
suboptimal control problem and Uncertain Systems are given. 

Chapter 2 Some of the current techniques for performance preserving controller 
reduction have been discussed in this chapter. The (P, 7 )-admissible reduction 
technique proposed by Goddard and Glover [12] and modified by Wang et. al. 
[26] have been discussed in some details. 

Chapter 3 The main results of this thesis are presented in this chapter. 
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Chapter 4 In this chapter the controller reduction techniques proposed in Chapter 
3 has been tested on the experimental HIMAT aircraft and a flexible launch 
vehicle. The complete yu-synthesis design of the autopilot for the launch vehicle 
in also included. 

Chapter 5 This concluding chapter summarizes the contributions of the thesis and 
points out the directions of possible further work in this problem. 



Chapter 1 
Preliminaries 


1.1 Introduction 

In this chapter we briefly describe the mathematical tools and the relevant background 
needed in our analysis of the Controller Reduction problem. Section 1.2 covers the 
definition of the generic state space system while Section 1.3 defines the Hardy Spaces. 
The concepts of Internal Stability, Linear Fractional Transforms and Coprime and 
Spectral Factorizations are outlined in Section 1.4, 1.5 and 1.6 respectively. The 
Standard Hoq suboptimal control problem is treated in Section 1.7 while Section 1.8 
covers Uncertain Systems. 


1.2 State Space Systems 

In this dissertation we only consider plants which may be described by a set of lin- 
ear, time invariant, finite dimensional, differential equations. Such systems may be 
represented by the following state space realization 


X = Ax + Bu 
y = Cx 4- Du 


( 1 . 2 . 1 ) 
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where x{t) e M"' is called the state vector, u{t) € is the input and y{t) € 
is the output. State matrices A, B, C, D are real, constant and of the appropriate 
dimensions. If p = m = 1 then the system is called scalar or SISO (single input 
single output) otherwise it is called multivariable or MIMO (multiple input multiple 
output). 

If the Laplace transforms of u(t) and y{t) are U{s) and Y{s) respectively and 
assuming zero initial conditions, the transfer function matrix from input to output is 
defined as G{s) where 

y(s) = G{s)U{s) (1.2.2) 


From this definition we have 

G(s) = C{sl - A)-'^B + D = 

Tab] 

where is used as a shorthand notation for G{s) 

-CD. 

Of all possible realizations (A, B, G, D) of G, those realizations with the smallest 
state dimension are called minimal. This is equivalent to (A, B) being controllable 
and (C, A) being observable. Given that (A, B, C, D) is minimal, the poles of G are 
defined as being the eigenvalues of A, and are closely related to the behaviour of the 
system. The system G is said to be stable if i?e(Aj(A)) < 0, and minimum phase if 
Re(z) < 0 where z is a transmission zero of the system. 

Consider a stable transfer function G with realization {A, B,C, D). Then the 
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be verified by direct substitution that P and Q satisfy the following Lyapunov 


equations 


AP + PA* + BB* = 0 
A*Q + QA + C*C = 0 


(1.2.5) 


pm'thernrore,(A, S) is controllable if and only if P > 0, and {C, A) is observable 
if and only if Q > 0. 

The eigenvalues of the product PQ are invariant under similarity state transfor- 
mations and hence form a set of input-output invariants. 


Pefiniti®’^ 1.2.1. Using the notation introduced above, the Hankel singular values 
of G ar6 defined as 

ai{G) = 


The cTj’s reordered such that they are descending in magnitude, ai > > cr„. 

The largest singular value gives the worst case or maximum gain of the system 
denoted by ^5 is widely used in Robust Control Literature. Detailed discussions 
on Linear Systems can be found in [2] and [4], while more on Hankel singular values 

can be found in [28]. 


^ 3 The Hardy Spaces 


For defining the Hoo space we require the following Hilbert Spaces. A detailed treat- 
ment of the Hardy Spaces can be found in [6] and [28], while Hilbert Spaces have 

been covered in [13]. 

Definition 1.3.1. £ 2 ( 1 ^) or simply £2 is a Hilbert Space of matrix valued(or scalar- 
valued) matrix functions F on jR such that 

/ H-oo 

trace[F*{jij)F{juj)]doj < oo 

-00 
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Definition 1.3.2. is a (closed) subspaee of £ 2 (jK)with the matrix functions F{s) 
analytic in Re{s) > O.The corresponding norm is defined as 
2 1 

IIFII 2 = sup — / trace[F*{a + ju))F{a + j(jS)\<ko 

a>0 27r 

Now we come to that class of functions bounded on the imaginary axis. 

Definition 1.3.3. /loo(jR) or simply Coo is a Banach Space of matrix valued(or 
scalar-valued) matrix functions F that are (essentially) bounded on jM,with norm 

Halloo = esssupa[Fijuj)] 

The rational subspace of Coo, denoted TZCoo consists of all proper rational transfer 
matrices with no poles on the imaginary axis. 

Definition 1.3.4. Tioo is a (closed) subspace of Coo with functions that are analytic 
and bounded in Re{s) > 0 i.e. in the open right half plane with the norm defined as 

ll-PIloo = sup a[F{s)] = suipd-[F{juj)] 

Re{s)>0 weR 

The Hoo space along with this norm is used to characterize transfer functions in 
Robust Control Theory. The rational subspace of Ttoo denoted by FfHoo consists of 
all proper and real rational stable transfer matrices. 

1.4 Internal Stability 

In the general feedback arrangement shown in figure 1.1, G and K are state space 
systems with the following equations 

w{t) 
u{t) 


x{t) = Ax{t) + [Bi B 2 ] 


(1.4.1) 
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Figure 1.1: The Generic Feedback System 


z{t) 

y{i)_ 


'Gi‘ 

G2_ 

x{t) + 

'Du 

D 21 

Di2 

D 22 

w{t) 

u{t)^ 


and K being described by 


(1.4.2) 


XRit) = AKXK{t) + (1.4.3) 

u{t) = CRXKit) + DKvit). (1-4.4) 

The system in Figure 1.1 is said to be well-posed if unique solutions exist for x{t), 
XK{t)i y{t) and u{t) , for all initial conditions x{0),xr(0), and all input functions 
w{t). Further such unique solutions must exist if the realizations are perturbed in 
some neighborhood of {A,B,C,D) and {Ak-,Br,Ck,Dr). 

Theorem 1.4.1 (Proposition 5.1 in [6]). The connection of G and K in Figure 
1.1 is well-posed if and only if I — is nonsingular. 

However well-posedness does not guarantee stability. 

Definition 1.4.1. The system in Figure 1.1 is internally stable if it is well-posed, 
and for every initial condition a;(0) of G, and (0) of K, the limits x{t), XK(t) 


hold when u; = 0. 
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Theorem 1.4.2 (Proposition 5.3 in [6]). The system of Figure LI is internally 
stable if and only if I — D 22 DK is invertible and 


-1 


^cl 

'a o' 

+ 

I 

-Dk ' 


' 0 

Ck 


0 


- —D 22 

I . 


. C 2 

0 . 


is Hurwitz. 


1.5 Linear Fractional Transforms 


This section defines the linear fractional transforms which shall be extensively used 
in this dissertation. See [28], [8], [27], and [6] for comprehensive treatments. 

Definition 1.5.1. Let M be a complex matrix partitioned as 


M = 


£{Pi+P2)'<{qi+q2) 


Mil ^12 

M21 M22 J 

and let A; e C<i 2 ^P 2 ^ other complex matrices. Then we define 

a lower LFT with respect to A/ as the map 


Ti{M, •) : 


with 

Ai) - Mil + Mi^Aiil - M22Ai)-^M2i 

provided that the inverse (/ — M 22 ^i)~^ exists. We can also define an upper LFT 
with respect to Au as 


with 

J^uiM,Au) := M22 + M2iAu{I — MiiAu) ^Mi2 
provided that the inverse (J — MiiAu)~^ exists. 
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A useful interpretation of an LFT is that !Fu{M, A) has a nominal mapping Mu 
and is perturbed by A, while M 12 , M 21 and M 22 reflect a prior knowledge as to how 
the perturbation affects the nominal map, Mu- This is why LFT is particularly useful 
in the study of perturbations, which is the focus of this dissertation. 

1.6 Factorizations: Coprime and Spectral 

In this section we briefly discuss the two type of factorization over TiTioo which will 
be the key tools in our new method of controller reduction. Details of both Coprime 
and Spectral factorizations have been discussed in [7], [28] and [6]. 

1.6.1 Coprime Factorizations 

Definition 1.6.1. Two matrices M and N belonging to “RHoq are right coprime over 
TZHoo ff they have the same number of columns and if there exists matrices Xr and 
1^ in TZHoo such that 

Similarly if two matrices M and N in TZTioo are left coprime over 'RTiao if they have 
same number of rows and if there exists matrices Xi and Yi in TlHoo such that 

[M N] = MXi + NYi = I. 

These two equations are often called Bezout Identities. 

Now let P be a proper real rational matrix. A right coprime factorization (ref) 
of P is a factorization P = NM~^, where N and M are right coprime over 72.77oo- 
Similarly we can get a left coprime factorization of P. 


= XrM + YrN = /. 
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Coprime factorizations are used extensively in both Ttoo synthesis and model ap- 
proximation. A particular advantage of such multiplicative decompositions for model 
approximation is the stability of the two coprime factors. 

1.6.2 Spectral Factorization 

We require the following definition. 


Definition 1.6.2. The conjugate transpose of G{s) is denoted 




■ -A* 

C* ' 

. -B* 

D* . 


Spectral factorizations may be viewed as the transfer function matrix equivalent 
of the square root. Given an Hermitian transfer matrix $ = which has no poles 
on the ju axis, spectral factorization allows us to calculate a stable, minimum phase 
spectral factor W such that $ = W^W . 


Theorem 1.6.1 (Chapter 7, Theorem 1, [7]). Every square matrix G{s) having 
the properties 

1. G,G ^ € 'K-Loo 

2. G-= G 


3. G{po) > 0 


has a spectral factorization of the form G = W'"W such that W,W ^ e TlHoo- 
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A matrix satisfying the conditions of the above theorem also has a co-spectral 
factorization G = WW^ such that W, W~^ G 'R^'Hoo 

1.6.3 Inner Functions 

Inner transfer functions are used extensively in both Hex, synthesis and model ap- 
proximation. 

Definition 1.6.3. A stable, p x m transfer function G, where p > m, is said to be 
inner if G~G = /. 

Definition 1.6.4. A stable, p x m transfer function G, where p > m, is said to be 
CO- inner if = I. 

The importance of inner functions stems from their norm preserving properties. If 
the input to an inner function has bounded energy then the output will be bounded 
by the same energy level. 

1.7 The Hoo Suboptimal Control Problem 

In this dissertation we aim to design reduced order controllers from full order con- 
trollers such that the closed loop structured singular value remains acceptable. We 
use the so called /i-synthesis along with the D — K Iteration to design the full order 
controller. The solution of the Tdoo suboptimal control problem, used as a part of the 
D — K iteration is presented below. Consider the general feedback structure of figure 
1.1. We use P here to denote the plant G to keep uniformity with usual nomenclature. 
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Here F(s) given by 



A 

Bi B2 

F(i>) = 

Cl 

Dll Di 2 


C2 

D21 D22 


where Bi 6 B2 € C\ € C2 € is a minimal state space 

realization of the generalized plant. The partition is assumed to be compatible with 
the signal dimensions. P{s) includes both the actual plant to be controlled and any 
frequency dependent weights used in the controller synthesis process. It is assumed to 
be FDLTI.The two Riccati equation solution was proposed in [ 5 ], and a most extensive 
treatment can be found in [ 28 ]. For LMI solutions one can refer to [6]. Earlier results 
can be found in [ 7 ]. 

1.7.1 Problem Statement 

Let 7 be some pre-specified performance level. 

Definition 1 . 7 . 1 . A controller K is said to be (P, 7)-admissible if K stabilizes P and 
WTzwWoo = m{p,K)\\^<j. 

The Hoo suboptimal control problem is to find a K that is (P, 7)-admissible. 
However, the solution also provides a parameterization for all controllers that are 
(P, 7)-admissible. 

The following assumptions are made for solving the problem: 

1. (A, B2, C2,D22) is stabilisable and detectable. 


2 . D22 = 0 
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0 

I 

1 D 21 = [0 /] and Dll ~ 

Dun 

Dun 

3. Dn = 




- Dim 

Din2 . 


where Dun € ]R{pi-"^2)x("^i-P2), Dun G £>^^21 G R^2x(mi-p2) 

and D1122 G R"^2xpa 


4. 


5. 


A - joji 

B 2 ■ 

Cl 

Dn . 

A - juil 

Bi ■ 

C 2 

^>21 . 


has a full column rank for all uj. 


has a full row rank for all u. 


1.7.2 The Suboptimal Solution 


Before stating the appropriate theorem it is convenient to define the following matri- 
ces: 

R = DlDi - Q j where Di = [Dn Du] 

and 

R = DiDl- where Di = 

Also assuming solution exists, define and Y^o as stabilizing solutions to the fol- 
lowing Algebraic Riccati Equations, 


Du 

Dn 


(A - BR-^DlCi)*Xoo + Xoo(A - BR-^DlCi) 

- XooBRr^B*Xoo + Ct{I - DiR-^Dl)Ci = 0 


(1.7.2) 


(A - BiDlR-^C)Yoo + Yoo{A - BiDlR-^C)* 

-YooC*R-^CYoo + Bi{I - DlR-^Di)Bi = 0 


(1.7.3) 
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and F and H are state feedback and output-injection matrices respectively, 

= ~R-^[DlCi + 

H = [HnHi2H2] 

= -[BiDl+YooC*]R-\ 

We are now in a position to state necessary and sufficient conditions for the existence 
of {P, 7 )-admissible controllers and a parameterization of all such controllers if any 
exists. 

Theorem 1.7.1 (Theorem 1 in [12]). For a plant described by equation 1.7.1 and 
satisfying assumptions 1-5: 

1. There exists an internally stabilizing controller K such that \\Fi{P,K)\\^ < 7 
if and only if, 

(i) ^ > max(a[DiniDiii2],a[Dl^-^^Dli2i]) 

(ii) there exists Xqq > 0 and Yoo > 0 satisfying equations (1.7.2) and (1.7.3) 
respectively such that pi^XooYaf) < 7 ^ 

2. Given that the conditions of part (1) are satisfied, then all rational internally 
stabilizing controllers K satisfying \\Fi{P, K)\\oo < 7 given by the lower LFT 
K — !Fi{Ka, ^),where $ € TZTtoo, such that 11^*1100 < 7 and Ka is given by 





Bal 

Ba2 


Kan 

Kal2 





Ka = 



~ 

Cal 

Ball 

Bal2 


. Ka2l 

Ka22 . 







Ca2 

Ba21 

0 


Ball — —Dn 2 iDliii{'y^I — DuiiDliii) - £>1122 




17 


where Dan S E"^ 2 xm 2 g are any matrices satisfying 

DanDl,2 = /-i5im(7"/-i5uii.-Diiiir'A*m 
DlnDan = / - ^>^112(7"/ - Aiiii 5 mi)~'l?iii 2 

and 

Ba2 = (-62 + Hi2)Dal2 
Ga2 = —Da2l{C2 + Fn)Z 
Bal — —H 2 + 5a2-Dai2-Dall 
Cal = F 2 Z DallDg^2iCa2 

I 

Aa = A + HC + Ba2D-,\Cal 

where 

Remark 1.7.1. i?i 2 a,-C* 2 ia, and hence i^a, are invertible. 

Remark 1.7.2. Aa—S 2 ajDi 2 a^ia are both stable, hence G 

1.8 Uncertain Systems 

In the last section we have given the primary result needed for synthesis of a subopti- 
nial controller for a feedback system where the plant model was completely specified. 
Now we seek to address uncertainty in the plant arising out of various reason ranging 
from incomplete knowledge of the parameters to dehberate under modelling. 
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The idea is to parameterize the entire set of uncertain plants using perturbations 
about a nominal mapping. 

We first present a motivating result from operator theory known widely as the 
Small Gain Theorem. 

Theorem 1.8.1 (Small Gain Theorem (Theorem 3.7 in [6])). Suppose Q is a 
member of the Banach algebra B. //'||(3|| < 1, then (/ — Q)~^ exists. Furthermore 

OO 

(/-<?)■’ = E«‘ 

fc =0 

1.8.1 The Structured Singular Value 

We define a commonly encountered structure in our uncertainty class; Let S and F be 
nonnegative integers, not all zero, and let n, si, ..., ss, /i, ..., Jf be positive integers, 
such that n = X) Si + X) /■!• Consider the subspace of n x n complex matrices 

C^s,F = {d^a5(<5iJsi, ...., <5s'Jss, Ai, ...., A^) : Si gC and Aj e 


and denote its unit ball by As^f = {A 6 CAs,f ■ <^(A) < 1}. Now define the time 
invariant uncertainty set as 

Ati = {A € C{C 2 ) : A is LTI and A(s) € As,f> for every Re(s) > 0} 

Definition 1.8.1. Given a matrix M € we define the structured singular value 

of M with respect to Ag^F by 


/i(M, As^f) = 


_1 

min{a-(A) : A 6 CA.s,f ond I — MA is singular} 


when the minimum is defined. Otherwise /r(M, As,f) is defined to be zero. 
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Computing n has turned out to be extremely difficult and efficient algorithms are 
still the subject of intense research. We have to be satisfied with an upper bound. 
For that we define the set of matrices © which commute with the perturbations, that 
is, ©A = A© for all A € ^s,F- Define the set of positive matrices 

^©s,F == {diag{Qi , ..., ©5, didfi, Ofifp) ; ©j € ©i > 0 and 6 i > 0} 
where H denotes the set of Hermetian matrices. 


Theorem 1.8.2 (Proposition 8.23 in [6]). 


IJ,{M, As,f ) inf a{Q^MQ 2) 

’ Q£VQs,f 


It can be shown that: 


Ati) = sup ia{M{juj), As^f) 

wGR 

< sup inf a(©^M©~2) 
wgR ©wSFSs.f 


1.8.2 Feedback Control of Uncertain Systems 


( 1 . 8 . 1 ) 


The generic feedback system is reprinted in Figure 1.2 but the perturbation A shown 
above incorporates the uncertainty in the plant while K represents the controller. The 
objective is to get such a K that will make the interconnection G and K internally 
stable which is equivalent to nominal stability of the interconnection. We also require 
that K should make the mapping 

= K), A) = Tu{M, A) 


small in some sense. 


Let M = 


Mil 
- M21 


Mi2 
M22 . 


be partitioned compatibly with the 


input and output signals. 
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Definition 1.8.2. With reference to figure 1.3 and with the extra assumption that 
the uncertainty class A is causal, the uncertain system (Mn, A) is robustly stable if 
(/ — MiiA)“^ exists in C{C 2 ) and is causal for each A € A. . 

Definition 1.8.3. With reference to figure 1.3 and with the extra assumption that 
the uncertainty class A is causal, the uncertain system (M, A) has robust performance 
if (Mil, A) is robustly stable and \\lFy,{M,A)\\ < 1, for every A € A. 

The next theorem gives the major result of ^-synthesis providing us with a test 
of the fj, norm to guarantee robust performance and stability. 

Theorem 1.8.3 (Proposition 9.10 in [6]). Suppose the configuration of Figure 
l.S is nominally stable, and the uncertainty class is Ay/ defined above. Define the 
perturbation set 

rr At, 0 1 

Ar/,p = \ : At, € Ay/, Ap € D,{C, 2 ) LTI, causal, ||Ap|| < 1 

I L 0 Ap J 

Then the following are equivalent: 

1. The uncertain system {M, Ati) satisfies robust performance: it is robustly sta- 
ble and ||.F„(M, At,)|| < 1, for every At, € Ati- 

2. The uncertain system (M, Api^p) is robustly stable. 

5- sup,^eKj^(M(iw), A5,ir+i) < 1 

We now shift our attention to the robust synthesis question. Prom the above 
discussion it can be shown that the robust synthesis for our standard setup in Figure 
1.3 can be recast as the finding of the infimum of 
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where K ranges over the stabilizing class 1C, and 0 ranges over the set of rational 
scalings Qti defined as 

= {0 e £(£ 2 ) : & nonsingular, LTI and 0(s) = diag{9i{s)I, ■..,9d{s)I)} 

The D - K iteration is a widely used heuristic to perform the above minimiza- 
tion.The steps are as follows: 

1 . Set 01 = I, 770 = 00 and the counter k=l. 

2. Solve for Kk in the Hoo synthesis infi^r^ '9k denote the 

achieved norm. 

3. Compare Vk with 7 ]fe+i;if they are approximately equal stop, and set K = Kk 
as the final controller. Otherwise continue. 

4 . Solve for 0fc+i in the scaled gain problem infe fc+1 
crement k and go to step 2. 



Chapter 2 


Controller Reduction: Performance 
Preserving Approaches 

2.1 Introduction 

The most widely used controller design techniques like Tioo loop shaping and 
synthesis typically produce very high order controllers. But practical controllers must 
be simple, linear and of low order. High complexity controllers are not only difficult to 
understand but they may turn out impossible to implement in hardware and software. 
Also from the point of integrity and reliability, simple controllers are preferable. 

The problem of reducing complexity via reducing the state dimension of the con- 
troller has been the subject of extensive research for the last forty years. See [1] 
and [11] for a complete list of references. In this dissertation we treat the controller 
reduction problem in a jU-S5Tithesis framework and derive a set of sufficient conditions 
for making the closed loop /r with the reduced order controller less than one. 

The design of reduced order controllers for high order plants may be divided into 
three categories [1] : direct design via optimization of controller parameters; model re- 
duction followed by the design of a low order controller using the reduced order model; 
and design of a high order controller from the high order model, followed by reduction 
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Figure 2.1: Different Paths to Controller Reduction 

of the high order controller. These three paths are shown in Figure 2.1. Each of the 
three paths may be split into two categories: firstly those dealing with strictly open 
loop considerations - making no reference to how stability or other appropriate closed 
loop objectives may be affected by the use of such reduced order approximations; and 
secondly those techniques which look at controller approximation in a closed loop 
framework. 

In this dissertation we will follow the third path and our goal will be to derive 
conditions that will guarantee closed loop stability and performance of the uncertain 
system with the reduced order controller. The approach will be to design a high 
order controller from a high order plant, followed by reduction of the high order 
controller. Information about the closed loop is incorporated into frequency weights 
and corresponds to requiring the modelling error to be smaller in selected frequency 
bands, e.g. around crossover. Hence the salient features of the high order controller 
are retained at the frequencies of most interest. 

Several Open Loop Approximation techniques form the basis of recent advances 
in the controller reduction literature. Balanced Truncation, a widely used technique. 
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is based on the elimination of the least controllable and observable modes. The 
balancing of the gramians, proposed by Moore[23] has been used by Glover[9] among 
others to find an upper bound for the approximation error in terms of Hankel Singular 
Values. 

The inability of this technique to provide any sort of optimal solution is over- 
come in the Optimal Hankel Norm Approximation Technique[9] which also provides 
a tighter upper bound in the approximation error than the balanced reduction. 

However there is one common problem for these reduction schemes, that is, the 
model to be approximated is required to be stable. The coprime factorization ap- 
proaches developed in [1],[20],[22] and [11], [12] have overcome this difficulty by de- 
composing the plant into two stable coprime factors, and then reducing the coprime 
factors instead of the plant. 

As has been discussed in great details in [1],[11] and [12], controller reduction is a 
distinct problem from model reduction. But since standard solutions are available for 
the open loop approximation problems, it is of practical convenience to express the 
controller reduction problem as a model approximation problem. This is achieved in 
this dissertation through frequency weights. 

This chapter is organized as follows: Section 2.1 discusses the three different paths 
to design a low order controller and mentions few of the open loop methods available 
in the literature. The frequency weighted Hankel norm approximation is covered 
in section 2.2 while sufficient conditions for stability preserving approximations are 
derived in section 2.3. The state of the art Tioo performance preserving controller 
reduction techniques are detailed in section 2.4. 



2.2 Frequency Weighted Approximation 

The problem can be defined as: 
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min 

K ^TZ.'Hoo 
deg{k)<d&g{k) 


Wi{K-k)W2 


oo 


( 2 . 2 . 1 ) 


where K € THi-oo, Wi and W 2 are the weights chosen to preserve the salient features 
of the controller at the frequencies of interest. One common method is the Weighted 
Balanced Truncation [28], where the balancing can be done from the solution of two 
Lyapunov equations. But the drawback of this method is that there is no a priori 
bound on the achievable frequency weighted error. So we intend to use the frequency 
weighted optimal Hankel Norm Approximation which is presented next. 


2,2.1 Frequency Weighted Optimal Hankel Norm Approxi- 
mation 


Since the design techniques to be presented in Chapter 3 lead to a two sided frequency 
weighted approximation problem with invertible weights, it is this framework which 
is now discussed. 

To calculate a solution via the Nehari Extension Theorem[7],[28] we solve the 
equivalent problem 

inf \\Mr\K-Q)Mr^^ (2.2.2) 

where Mi, M 2 , € TZHoo such that M 2 — W 2 W 2 and MiM^ = 
W^Wi, and show that such a solution satisfies bounds on equation 2.2.1. 

Define H = [M^”'^i<'M^”^]+ E TZTCoo- Then calculate an optimal Hankel norm 
approximation H of H, 


inf 

HeTlH] 


00, (k) 


H-H 


H-H 


H 




(2.2.3) 
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To calculate the required approximation to K notice that 
Q = M2~(H + e 

Therefore setting K = [file, the stable part of Q gives a stable )=''> order frequency 

weighted approximation to K . 

The following lower bound or error for frequency weighted approximation can be 


derived: 


< _inf \W 2 {K-K)Wi 

deg{k)<deg{K) 


(2.2.4) 


where e IbWos has McMillan degree „ and X £ KtU haa McMUlan degree 
k(k < n). Ideally we would like reduction techniques with an upper bound close to 

this lower bound. 

An upper bound on the error for weighted Hankel norm approximation is much 
more difficult to obtain than for the unweighted case. Anderson[l] derives a bound 
for weighted Hankel norm approximations which is proportional to the unweighted 
error. The value of the factor of proportionality is dependent solely on the weighting 
functions. Glover et al[10] present the following error bound, 


W2{K - K)Wi 


<ak+i{H)+ \[H]--W2 


However, if either of the weights have zeros close to the imaginary axis then the second 
term can become arbitrarily large. Since this bound cannot be calculated a priori, 
(nor is it generally tight,) it is not as useful as the error bound in the unweighted 
case. For practical purposes, a combined model reduction and parameter optimization 


approach is used. 
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Figure 2.2: The generic Feedback System 


Now we turn our uttention to derive weights leading to open loop frequency 
weighted approximations, such that certain closed loop properties ate retained. For 
a start we present the sufficient conditions for the weights to be stabilizmg. 

2.3 Stability Preserving Weights 

We again consider the generic feedback system redrawn below in Figure 2 . 2 . where 
G is given by 

Gii Gi2 

G = 

. G21 G22 - 

and <?22 = ^ is a p X g transfer matrix. Suppose if is a order 

1 C 2 D22 J . t.- 

controller that stabilizes the closed loop system. We ate interested in investigating 
controller reduction methods that can preserve the closed-loop stability with ' the r^ 
duced order controllers. The following results are for the so called additive reduction. 

Lemma 2.3.1 (Lemma 19.1 in [28]). Let K be a siuMimng controller and K be a 
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reduced order controller. Suppose K and K have the same number of right half plane 
poles and define 

A:=k-K, Wa:={I-G22K)-^G22 
Then the closed-loop system with K is stable if either 

IIW^aAlloo < 1- (2.3.1) 

or 

IIAW^ailoo < 1. (2.3.2) 

In view of the above lemma, the controller K should be reduced in such a way 

so that the weighted error Wa{K — K) or {K — K)Wa is small and K and 

OO CO 

K have the same number of unstable poles. If K is unstable, K is usually separated 
into stable and unstable parts as 

K = K+ + K. 

where K+ is stable and only K+ is reduced to K+, and the final reduced order 
controller is given hy K = K+ + K-. 

2.4 7-^00 Performance Preserving Techniques 

In this section, we present a brief summary of current techniques od controller re- 
duction preserving Hco performance. The assumptions of Hoo-suboptimal solution 
(Section 1.7.1) apply also to this section. 

A major approach in this respect is the approximation in gap metric. However 
we choose to skip it here in order to concentrate on approaches which motivate our 
development in the fi framework. A bound on the performance degradation was 
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given by Lenz. et. aL[19]. When K = + W 2 AW 1 is a nominal (P,7) admissible 

controller subject to a weighted perturbation, they proposed a method to calculate 
weights such that 

\\W^\K - (2.4.1) 

Equation 2.4.1 represents a sufficient condition for controller K to stabilize the closed 
loop and give a performance bound of \\Ti{P, K)\\^ < \/2. This represents a (possi- 
ble) degradation in performance of y/2. A similar result may be derived by assuming 
K is & nominal controller subjected to coprime factor perturbations. 

This bound was improved and also a practical algorithm for calculating the weights 
was proposed by Goddard and Glover in [11], [12]. 

2.4.1 (P, 7 )- Admissible Controllers 

As with the Lenz approach[19], in this case the sufficient conditions are expressed 
in terms of weighted error between a (reduced order) controller K and a (nominal) 
(F, 7 )-admissible controller e.g. 

\\W^\K - K°)W{-^\^ < 1 ^ \\HP,K)\L < 7- (2.4.2) 

For the coprime factor perturbation approaches, the sufficient conditions axe expressed 
in terms of the weighted error between left coprime factors of a (reduced order) 
controller and a nominal (P, 7 )-admissible controller, e.g. 

l|lV2-H[t/ - IC/" 1 ^ Jfjlloo < 7 . (2-4,3) 

where K = V~^U and Equation 2.4.2 and 2.4.3 have a natural 

interpretation as frequency weighted model reduction problems. Given appropriate 
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weights, the goal is to calculate reduced order K or \U V], which lies within an Cqq 
ball (of right radius Wi and left radius W2), of or [U^ K is then a reduced 
order (P, 7) -admissible controller. The larger the £00 ball easier the approximation. 
We therefore wish to make Wi and W2 large in some sense. 

Assume that any controller K may be represented by a nominal admissible con- 
troller subjected to a weighted additive perturbation i.e. K — + W2AW1, where 

Wi and W2 are stable minimum phase and invertible and A is stable. This assump- 
tion imposes the restriction that K and possess the same right half plane poles, 
thus to some degree limiting the set of reduced order controllers attainable by this 
method. 

Consider Figure 2 . 3 . Observe that P/(P, K) = A) and hence 

miP, A) 1 U < 7 A) < 1 ( 2 . 4 . 4 ) 

00 









where 
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I ] I 7-57 0 

) 

j 0 J L 0 W 2 

<S(-,-) denotes the Redheffer star product. Suppose is some nominal (R, 7 ) ad- 
missible controller. If stable, minimum phase and invertible weights Wi and W 2 can 
be calculated such that R is inner and R 21 has full row rank (note that this is not 
generally possible-see Lenz et al [19] then 

liAIL < 1 < I 

^ <1 (2.4.6) 

00 

« ||/■,(p,^r)lu<7. 

Although it is not possible to calculate weights such that the above necessary and 
sufficient condition is satisfied, it is possible to calculate weights which make R a 
contraction. Using such weights gives the sufficient condition, 

\\W^\K - K^)W,-^\\^ < 1 il^K-P,^)lloo < 7- (2.4.7) 


(2.4.5) 



Using similar methods the sufficient condition on the coprime factors as in equation 
2.4.3 can be derived[ 12 ]. A method to calculate the weights so as to make R a 
contraction is given in the following theorem. 


Theorem 2.4.1 (Proposition 4.5.1 in [11]). Assume j|i?iil| 7 define 
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Define X = {W^Wi)~'^ and Y = (yV^W^Y^- Then R is a contraction if and only if 


■ A 

0 ■ 

> 

■ Li 

L2 ■ 

. 0 

y . 


- T 2 

Lz . 


To obtain the set of controllers K such that 1|A1 |qo < 1 is equivalent to finding 
the set of controllers K that lie within a weighted £oo ball of left radius W 2 and right 
radius Wi centered on K^. Hence we require both Wi and W 2 to be large in some 
sense to make any subsequent approximation as easy as possible. In the scalar case 
maximizing |H/ik 72 l is appropriate, however in the matrix case the optimal choice is 
not clear. Goddard and Glover [11], [12] have proposed maximizing the the product 
trace(Wf^Wi)trace(W 2 W 2 ), and the product det(Wf^Wi)det{W 2 W 2 ). We state only 
the relevant theorems giving the results of the respective optimizations. 

Lemma 2.4.2 (Lemma 4.6.1 in [11]). Define L as in Theorem 24 . 1 and assume 
Li > 0, with trace{Li) 4 0,and L 3 > 0. Also define 

aopt = y/{trace{Ls)jtrace{Li) 

Zl = L2L2 
Z 2 = L2Zfi^L2 

where Zi = Zf' > 0 and Z 2 = Z 2 >Q- Then the weights Wi and W 2 which maximize 
the product trcuce{Wi^Wi)trace{W 2 W^) subject to equation 24-9 satisfy 

{WjfWi)-^ = Li + a^tZ 2 (2.4.10) 


(W 2 W 2 ~)"^ = Lz + aoptZi 


(2.4.11) 
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Lemma 2.4.3 (Lemma 4.6.2 in [11])- Define L as in Theorem 2.4.1 and assume 
Li > 0 and Lz > 0. Also define SS'" = Li > 0, T~r = L 3 > 0 , 

Z 2 r r r n'^—l 

I — O JU 2 J^^ 1^2 ^ 

where Z\ — > 0 and Zi = > 0. Then the weights Wi and W 2 which maximize 

the product det{WiWi)det{W 2 W 2 ) subject to equation 2.4.9 satisfy 

{Wfi'Wi)-^ = S{I + = Li+ (2.4.12) 

(WaWT)'^ = r~(J + Z2)T = Lz + T^Z2T. (2.4.13) 

A slightly modified approach has been proposed by Wang et al in [26], which 
gives lower order weights than the previous technique. Instead of using additive 
perturbation on the controller, this method is based on additive perturbation on the 
closed-loop transfer function. 

Consider the class of closed-loop transfer functions with reduced order controller 
which can be represented as 

Gr = G- W 2 AW 1 

where G is the closed loop transfer functions with full order controller, A is stable 
perturbation, Wi and W 2 are stable, invertible minimiun phase weights. The following 
theorem gives the sufficient conditions for performance preserving reduction. 


Theorem 2.4.4 (Theorem III.l in [26]). Assume ||G||oo < 1, awd W 2 are 
stable, invertible and minimum phase weighting functions and also satisfy 


I 


W2W2 


G 


G 

> 0 

I-W^Wi . 


(2.4.14) 
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\\W^\G-Gr)W{^\\^<l - (2.4.15) 

Then 

ll*^r||oo < 1 

Remark 2.4.1. From the above theorem it can be seen that W^Wi and W 2 W 2 should 
satisfy the following inequalities: 


I - W2W£' > Gil - W^Wi)-^G 


I - W{^Wi > G^(I - W2W2T^G 

Since Kr is obtained by frequency weighted model reduction, it is preferable 
that and are as small possible in some sense. Instead of the maximiz- 
ing traceiW^\Vi)traceiW 2 W 2 ), Wang et al[26] chose to minimize 


trace[W2^GW{^W{^^G^W2^^] 

The following theorem gives the sufficient conditions. 


Theorem 2.4.5 (Theorem III.2 in [26]). < 1 and 


I - W 2 IF 7 G 

G- I- Wf W - 1 


> 0 


(2.4.16) 


then IIVF 2 ^GW^ reaches the minimum ortTace]yV 2 ^GW-^ reaches 

the minimum with W^Wi = / — y/G'^G and W 2 W 2 = I ~ VGG^ 


2.5 Conclusion 

In this chapter we have reviewed some of the current techniques for performance 
preserving controller reduction. We have left out all the open loop approaches except 
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Hankel norm Approximation in order to emphasize on topics motivating the next 
chapter. The (P, 7)-admissible reduction technique proposed by Goddard and Glover 
and modified by Wang et. al. are especially important as they motivate our derivation 
of the sufficient conditions for performance preserving reduction in the p. framework 
to be presented in the next chapter. 



Chapter 3 


Controller Order Reduction of 
Uncertain Systems 

3.1 Introduction 

In the last chapter we reviewed several controller reduction techniques which guar- 
anteed Hoo performance of the closed loop system with the reduced order controller. 
These methods in general do not consider any structure in the uncertainty. However 
consideration of structure in the uncertainty may make the reduction schemes less 
conservative. Thus a greater reduction in the controller order may be possible. More- 
over the proposed techniques integrate closely with the /x-synthesis controller design 
framework. 

In this chapter we present two new controller reduction schemes which preserve the 
closed loop jLi value to be less than one, thus guaranteeing closed loop robust stability 
and performance. Consideration of structure in the uncertainty is thus embedded in to 
the reduction algorithms and a generalization of current controller reduction schemes 
for the /X framework is achieved. Section 3.2 develops a new proof for Kavranoglu’s 
[17] additive reduction on the full order controller while 3.3 proposes a coprime factor 
reduction algorithm. 
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Figure 3.1: The Generic Feedback System with full order controller 

3.2 Additive Perturbation Reduction 

Let us again consider the generic uncertain feedback system of Figure 3.1. The 
uncertainty € Ayj as defined in Theorem 1.8.3. We use the same definitions of 
CAs^F, ^Tl,p as in Section 1.8. Let if be a controller that makes the closed 

loop uncertain system (!Fi{G,K),A) robustly stable and having robust performance 
according to Definitions 1.8.2 and 1.8.3. Thus by Theorem 1.8.3 

/i(M, Ati,p) = sup n{M{juj),As,F+i) < 1 (3.2.1) 

cjeR 

where M = Ti{G,K). Conversely we can say, again from Theorem 1.8.3 that for 
any controller K, if we can show that ii{M{juj),As^F+i) < 1 then the uncer- 
tain feedback system (M, Axi) is robustly stable and has robust performance. This 
observation is the key to to the following controller reduction procedure. 

Next we present a result that will be the basis of our controller reduction schemes. 
The main objective of the following result is to extend the standard mixed- jj, analysis 
to the case when the complex uncertainties have frequency dependent upper bounds. 
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We augment the uncertainty structures to meet the mixed-/^ setup. Let R, S and F be 
nonnegative integers, not all zero, and let n, ri, ..., r/j, si, ..., ss, /i, /f be positive 
integers, such that n = + X] + X] /i- The frequency dependent bounding set 

is thus defined as 

W = {W = diag{piIr„..;pRlrR,WiIsi,--u)shs,'^S+i^fi,-,'uJS+Flfp) 

: Pi > 0,i = 1,2, ...,R, and some functions (3.2.2) 

Wi ■ M — > 1^+) ^ ~ 1) (S' T -F)}- 

Let CAre be the subspace of real x Jfri matrices defined by 

CAre = {diagijilr^, ..., 'YRirn : 7i G M} (3.2.3) 

Next we define a structure consisting of those functions that have continuous extension 
on the right half s-piane. Define the set: 

CAxi = {A(s) € Hoo ■ A is LTI and A(s) € CAs^f-, for every Re{s) > 0} 

(3.2.4) 

Next we define the augmented structure 

Arc = {Arc : Arc € diag{Ar, A), Ar € CAre, A € CAtj} (3.2.5) 

Finally we define the uncertainty set bounded above by the frequency dependent 
functions. Given W € W, define 

Aw. = {A € Arc : A(ia;)*A(ia;) < W{ojf, 'iuj € M} 

Then we have the following result 

Lemma 3.2.1 (Theorem 6 in [25]). Let W eW and G € Ifoo- {f 

'fa(Wiuj)G{jcu),Aw) <1, Vo; € R, 


then (C, Ary) is robustly stable. 
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Figure 3.2: Uncertain system with the reduced order controller 


We are now ready to give our main result for performance preserving controller 
reduction. To reframe the performance preserving controller reduction problem in 
the fi setting, consider the following block diagrams. (Figure 3.2 and 3.3). Let the 
reduced order controller Kr be expressed as a perturbation to the full order controller 
Kf such that Kr = Kf + Ac. We assume that Ac € For the upper bound we 
define a frequency dependent function Wa = 7 (u’)J where 7 (w) : R — >• R+. Clearly 
Wa e W. Now define 

T = S{G, ^ ) (3.2.6) 

. I 0 . 

where <S(-, •) denotes the Redheffer star product. 

- r ^ 0 ■ , ^ 

T = T (3.2.7) 

.0 Wa . 

and 

Aupc = ( ■ ^Ti,p e ^Ti,p o-nd Ac € Aw i (3.2.8) 

I 0 Ac J 
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Figure 3.3: Plant and controller interconnection for additive reduction 


A main drawback of this method is that this characterization of the set of reduced 
order controller limits Kris] to have the same number of unstable poles as A/(s). 
This problem is addressed in the next section. 

Proposition 3.2.2. If a{Ac{ju})) < 7(w), Vw, and iJ,{Tijuj),Aupc) < 1, Vw, then 
the closed loop uncertain system itFiiG,Kr),Axi) is robustly stable and has robust 
performance. 

I 0 

Proof. Clearly G W and by hypothesis for each oj 

0 Wa 

a(Ac(jw)) < 7(0;) =» Ac(ja;)*Ac(jai) < Wf 

r 1 ^ 

I 0 

^ Aupcij^] ^UpciJ^) 

0 Wa 


I 0 

Now T = 

- 0 Wa 


T. Also by hypothesis we have ;u(T(jw), A„pc) < 1- Thus 
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by direct application of Lemma 3.2.1 we have the system (T^jcj), A^pc) is robustly 
stable. Now for any particular Kf and Kr we have 


{T(juj),Aupc) ={S{G{juj), 


Kf 

0 


0 

/ 



) 


= {Ti{G{jaj),Kr),ATi,p). 

Thus we have shown that the closed loop uncertain system {J^i{G, Kr), Axi,p) is 
robustly stable. But by theorem 1.8.3 this is equivalent to {J^i{G,Kr),Axi) being 
robustly stable and having robust performance. Hence proved. □ 


Remark 3.2.1. Prom the above theorem it is evident according to Theorem 1.8.3 
that if we can find an Wa such that ij,{T{juj), Ay,pc) < 1, Vca then any Kr satisfying 
d'(Kf — Kr) < 7(a)) is a robustly stabilizing controller and also provides robust 
performance to the closed loop system. 

Now let 7(ja)) be a rational function approximating 7(0;) such that j(ju>) < 7 (a>). 
Thus for each cu we have the following set of equivalences. 

a{Kf{ju) - Krijoj)) < 7(iw) 

^(ju))^iKf{jaj) — Kr{ju!)) < 1 
^ a[^^I{KfUu^) - KrUu^))] < 1 
« a[W-\Kf{ju:) - Kr{juj))] < 1 
^ \\[W-\Kfijio)-Kr{ju))]\\^<l 

Thus we have converted the controller reduction problem into the widely studied 
frequency- weighted Coo model approximation problem. This problem can be solved 
by a number of methods such as frequency weighted balanced truncation and optimal 
Hankel norm approximation. 
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Now the remaining problem is to calculate Wa such that jLi(f, Aupc) < 1- This is 
done by the following technique. We have already assumed that Wa{oj) — 7(w)/. We 
can find 7 (w) using among others the /j, — tools software of MATLAB [8]. One has to 
perform a few bisection steps on the size of a{Kf{ju}) — Kr{ju))) for each frequency 
to determine the maximum size possible such that fi{T{joj), Ay,pc) = 1. 

Next we show that when expressed as a two sided frequency weighted approxi- 
mation problem such as ||Wi(Ary — i<rr)W 2 ||^ < 1, the Wi and W 2 produced by the 
above algorithm minimizes the weights in the sense that trace{WiWi)trace{W 2 W£') 
is min im ized. For the next result we need the following standard lemma [26]. 


Lemiha 3.2.3. Suppose a Hermitian matrix is partitioned as 


A B 
B* C 


where A = A* 


and C = C** > 0. This matrix is positive-semidefinite if and only ifA—BC ^B* > 0. 


We have the following proposition. 

Proposition 3.2.4. If Wi and W 2 be two stable, invertible weighting function with 
minimum phase, satisfying Wf'Wi = and = l{j^)I then 

\\Wl{Kp - Kr)W2\\^ < 1 ^ HKfijcv) - KriM) < %ju>) Vcu 


and trace{W^Wi)trace{W 2 W 2 ) is minimized. 



44 


Proof. Let M(s) = Kf{s) — Kr(s). We have 


IIW1MW2II < 1 

^ WiMW2Wf'M'^Wf' < I 

MW2W2~M~ < (Wf Wi)-i 

(WfWi)-^ M 

M~ (W2W2~)-^ 

1 

0 M 
[M~ 0 

7/ M 




> 0 


M~ 7/ 


> 


> 


0 


-(W2W2)-^ 


0 

7 /-(WfW 3 )-^ 0 

0 7 / - (W 2 W 2 ~)- 


Now if WiWi > j(ju))I and W2W-^ > j{ju)I, we have 

7/-(WfWi)-i 0 

0 7/ - (W2W2~)-^ 

7I M 

> 0 

M~ 7/ 

M~M < 7^/ 

Therefore at each frequency 


>0 


o-(M(iw)) < j{juj) 

a{Kf{joj) - Krijoj)) < ^{juj) 'ibj 

Under this conditions trace{Wf'Wi)trace{W2W2) is trivially minimized for Wf'Wi = 
W2Wf' = f{juj)I. Hence proved. □ 
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3.3 Coprime Factor Reduction 

3.3.1 Introduction 

One of the main drawbacks of the above method is that it assumes the controller to 
be approximated is stable. Robust controller design techniques like Hoa loop shaping 
or yu-sjmthesis do not guarantee stability of the controllers. However it is possible to 
successfully approximate unstable systems. One simple method, if K is unstable, is to 
separate K into stable and unstable parts as K = K+ 4- K- where K+ is stable. Only 
K+ is reduced to K+, and the final reduced order controller is given hy K = K+-\-K-. 
If K has jui axis poles then they may also be copied directly into K but must have 
precisely the same residues. 

The implication here is that no approximation can be undertaken on the unstable 
part of the controller. The reduced order controller is therefore constrained to possess 
the same unstable poles as the original controller. However, there is no reason why 
good approximations should possess this property. 

An alternative to the above procedure is to consider a fractional description, the 
basis of which is multiplicative decomposition rather than additive decomposition. 
Fraction approaches to controller reduction have received much attention in the lit- 
erature. These are generally termed as coprime factor approaches. Various sufficient 
conditions for the stability of the closed loop system when the controller is subjected 
to coprime factor perturbation has been derived in [20], [22]. Tioo performance pre- 
serving techniques with multiplicative uncertainty have been studied in Goddard and 
Glover [12]. 
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3.3.2 Coprime Factor Approach 

In this section we propose a coprime factor based controller reduction scheme that 
keeps the closed loop /u value with the reduced order controller to remain less than 
one, thus guaranteeing robust stability of the closed loop system. Apart from having 
the natural advantages of an multiplicative scheme as mentioned m section 3.3.1, we 
also show that the perturbations have a block diagonal structure. This points to the 
possibility of greater reduction in the controller order than possible with the additive 
scheme described in Section 3.2. 

We assume that the left coprime factors of Kr = V-^Ur may be represented by 
perturbations to the left coprime factors of the full order controller Kf = Vf% i.e. 

Ur = Uf + Au (3.3.1) 

W = Vf + Av (3-3-2) 

We assume A^;, Let = jv{^)I and Wu = wW where 7y (a;), 7tr(^) •• 

R K+. We use the same definition of W as in Section 3.2. Clearly I^y , Wu € W. 
Figure 3.2 can be redrawn as in Figure 3.4, which after some manipulations give 

Figure 3.5. We define 

'Kf V 

KfA= I 0 0 (3-3-3) 

-Kf 

T = S{G,KfA) 


and 


/ 0 0 

f = 0 Wy 0 T (3-3-5) 

0 0 Wu 
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Figure 3.4: Uncertain system with reduced order controller 



Figure 3.5: Plant and controller interconnection for coprime factor reduction' 
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while 


^upvu ~ 


■ ^Ti,p € Aj'j.p, Ay € Aw, Au e Aw> (3.3.6) 


Ati,p 0 0 

0 Ay 0 
0 0 Ac; 

We have the following proposition to motivate our controller reduction procedure. 

Proposition 3.3.1. Let d-(Ay(ja;)) < yviuS) and a{Ajj{juj)) < 7c;(w), Vw. Now if 
fj-{T(juj),Aupvu) < 1) Vw, then the closed loop uncertain system {fFi{G,Kr),A'pi) is 
robustly stable and has robust performance. 


Proof. Clearly 


J 0 0 

0 Wy 0 
0 0 Wc; 


€ W and by hypothesis for each w 


a-{Av{joj)) < 7v(w) and a{Au{joj)) < 7 c;(w) 



1 

o 

* 

Ay 

0 


Wv 

0 

4=^ 

0 Ac; 


0 

Ac; 

< 

0 

Wu 


Aupvuij^') < 


J 0 0 

0 Wy 0 
0 0 Wc/ 


(3.3.7) 


Now T = 


/ 0 0 
0 Wy 0 


0 0 Wu 

by direct application of 


T. Also by hypothesis we have p{T{ju), Aupvu) < 1- Thus 
Lemma 3.2.1 we have the system (T(jc<;), Aypt,„) is robustly 



stable. Now for any particular Kj and Kr we have 


= 5 G(ia;), 


Kf 

I 




-Kf -V, 


-1 


Vf^ 


■‘T’ 


I ^upvu 


(3.3.8) 


= {:Fi{G{ju),Kr),ATi,p). 

Thus we have shown that the closed loop uncertain system (J^i{G,Kr),ATi,p) is 
robustly stable. But by theorem 1.8.3 this is equivalent to (J^i(G, Kr), Atj) being 
robustly stable and having robust performance. Hence proved. □ 


Like in the additive case the above result implies that if we can find suitable 
7 v(w) and such that Aupvu) < 1 ) then our controller reduction 

problem can be written as frequency weighted Too approximation problems. Let 
and 7 i 7 (jw) be two rational functions such that Wv{joj) = < 7 k(^) 

and Wu{ji^) = < 7 i 7 (^t^)Vw. Then our controller reduction problem reduces 

to the following set of frequency weighted Too approximation problem following the 
arguments outlined in the additive case. 

\\Wy\V^-Vr)\\^ < 1 (3.3.9) 

\\W^HUf-Ur)\\^ < 1 (3.3.10) 

Thus our objective again reduces to finding 'yv(ui) and jui'^) such that for each fre- 
quency Aupvu) < 1- Here also the bisection algorithm outlined in the addi- 

tive case may be used. But search over two separate parameters is not practically pos- 
sible as the dep enhance of fi on the bound of the uncertainty is not explicitly known. 
However we could just search over one parameter by taking Wy = Wjj — jyu{uji)I. 
Under this assumption we find out 71 / 17 ( 0 ;) by performing a few bisection steps on 
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the size of a 


Av(iw) 0 


for each frequency to determine the maximum 


0 

size possible such that Aupvu) = 1 Vw. The weighted Too approximation 

problem can be solved by a number of available algorithms such as weighted bal- 
anced truncation or optimal Hankel norm approximation. In omr design problems 
discussed in the next chapter we use the frequency weighted optimal Hankel norm 
approximation for its tighter upper bound of approximation error. 


3.4 Summary 


In this chapter two performance preserving techniques in the //-framework has been 
proposed. A new proof has been given for the Kavranoglu’s additive reduction tech- 
nique and the weights for two sided reduction have been derived. A new coprime 
factor reduction scheme has been proposed which guarantees closed loop stability 
and performance with structured uncertainty. The coprime factor perturbations to 
the controller have been found to have a block diagonal structure thus improving the 
reduction algorithm. 



Chapter 4 

Design Examples 


4.1 Introduction 

In this chapter we deal with the practical design aspects of /i-synthesis and controller 
order reduction techniques discussed in Chapter 3. The two numerical examples 
consider ^-synthesis of reduced order controllers for LTI systems. Given a generalized 
plant P and the uncertainty set A , the objective of //-synthesis design is to find a 
controller K which stabilizes the closed loop uncertain system Pu{P, A) and makes 
the norm WPuiP^iiP, ^)li < 1- The generalized plant incorporates both the actual 
plant which is to be controlled and weighting functions which specify performance 
objectives. The performance weights are selected to emphasize objectives such as 
disturbance attenuation at low frequency and robustness to modelling uncertainty. 
Selection of appropriate weights is nontrivial and often involves a certain amount of 
iteration. For a more detailed discussion on weight selection see [8], [3] and [18]. 

In Section 4.2 we examine the application of the proposed algorithms on the 
pitch control of an experimental aircraft model. This example was chosen as it has 
been widely studied' in the controller reduction hterature and thus can be used for 
comparative purposes. Section 4.3 details the complete (i synthesis design and the 

." T-.Trj 

iro A 
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d4 

Figure 4.1: The HIMAT Closed Loop Interconnection 

reduction of the full order controller for the experimental data available from a flexible 
Launch Vehicle. 

4.2 Pitch Control of the HIMAT Model 

4.2.1 The Plant Description 

In this example the controller order reduction algorithms of Chapter 3 are applied to 
the design of a pitch axis controller of an experimental highly maneuverable airplane, 
namely HIMAT. A linearized model and detailed specifications can be found in [8]. 

A diagram for the closed loop system which includes the feedback structure of the 
plant and the controller, and elements associated with the uncertainty and perfor- 
mance specifications, is shown in figure 4.1. The actual aeroplane is represented by 









53 


the dashed box with the transfer function G, with the combination of and Aq 
representing the amount and type of uncertainty in the nominal model Gnom- Gnom 
has four states 

5u-perturbation along the velocity vector 
a-angle of attack 
g-pitch rate 
0-pitch angle 

two control inputs, 

0e-elevon command 
^c-canard command 

two measured outputs, 

a-angle of attack 
0-pitch angle 

The flight path angle is given by p = 0 — a. In [8], a 30 state Hcc controller is designed 
using p-synthesis to give robust performance for the following three maneuvers: 

Vertical Translation: Control the vertical velocity at constant 0 (ct varies). This 
implies that the attitude remains constant as the velocity vector rotates. 

Pitch Pointing: Control the attitude at a constant flight path angle (i.e. 6 — a). 
For this maneuver the velocity vector does not rotate. 

Direct Lift: Control the flight path angle while holding the angle-of-attack con- 
stant (i.e. p = 0). This maneuver produces a normal acceleration without changing 
the angle-of-attack. 
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d2 
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d4 


Figure 4.2; The Augmented Plant and Controller 


These performance objectives are captured in the weight Wp. Thus the design 
objectives are translated into the requirement that the transfer functions from the 
disturbances [di ds d^]'^ to [ei 62 ]^ to be less than unity for all i| Ag|1 < 1. 

A generahzed plant incorporating all the weights is first formed. The resulting 
plant Gaug along with the uncertainty set Ag and the controller Kf is shown in 
figure 4.2. 


4.2.2 Controller Reduction Results 


The generalized plant Gaug h^s 10 states. AT/ is a 30 state controller satisfying 
the performance specifications. The closed loop robust performance is conveniently 
measured in terms of the p of the closed loop with respect to an augmented uncertainty 
set. We define 


Ai 0 
0 A 2 


: Ai ec2^^A2eC' 


'2x2 


^AUG ■= 
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Reduction Approach 

Controller Order 

fiiPi{P,K),A) 

Full Order 

30 

0.966 

p^A{P) 

16 

1.0567 

mC(P) 

15 

1.0536 

M(P)+b.t. 

8 

1.1219 


Table 4.1; Results of Model Reduction 


According to theorem 1.8.3 the Augmented plant Gaug has robust performance and 
robust stability if iJ^{J^i{GAUG,Kf),AAUG) < 1- In this case Kf is found to satisfy 
IJ,{J^i{GAUG,Kf),AAUG) <0.996 

The results of model reduction are presented in Table 4.1. The following abbrevi- 
ations are used 

IJ.A{P) Additive Perturbation Reduction 
liC{P) Coprime Factor Perturbation Reduction 

This can be compared with the results from [12] and [26]. Table II from [26] is 
reprinted below. In Table 4.2, A{^), CP{^) and CP{e) are the methods proposed 
in [12]. UOHNA represents unweighted optimal Hankel norm Approximation, UBT 
represents unweighted balanced truncation and LCFR represents unweighted optimal 
Hankel norm reduction of controller left coprime factors. 

Thus the proposed algorithms are found to reduce the 30*^ order controller to 
comparable orders as in the table above. It should be noted that the full order 
controller used in Table 4.2 above is of order 20 while the proposed approaches in 

Table 4.1 uses a 30*'^ order full controller. 

It is suspected that further reduction is possible but for the poor performance of 
the actual frequency weighted reduction algorithms. In fact a 8*^ order controller 
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Reduction Approach 

Controller Order 

mp,K)\L 

Full Order 

20 

0.97 

A(4) 

14 

0.99 

CP{^) 

16 

0.98 

CPiO) 

16 

0.99 

UOHNA 

13 

0.98 

UBT 

13 

0.98 

LCFR 

13 

0.98 

A($) 

7 

1.27 

CP{^) 

7 

7.07 

cpIq) 

10 

2.02 


Table 4.2: Results from [26] 


can be found simply by balanced truncation of the reduced controller of the ij.A{P) 
approach, which satisfies the sufficient condition i.e. ii{Ti{GAUG,Kr).,/^AUG) < 1- 
It should further be noted that the main advantage of the proposed methods is 
not utilized in this example as the original uncertainty of the HIMAT system was 
unstructured. With the incorporation of structure in the uncertainty the proposed 
methods remain applicable. This is not the case with most of the currently available 
techniques which do not incorporate any knowledge of the uncertainty structure. 

Figures 4.3, 4.4 and 4.5 show the closed loop ^ plots for the full order and the 
reduced order controllers obtained by each of the proposed algorithms. The peak 
value Of 0.966 for the full order controller increases to 1.05 for the reduced order 
controllers. However for further reduction the peak [j, value increases to unacceptable 
high values. So 16 with M(^) with ij,C{P) are the minimum orders achievable 

with the proposed reduction techniques. 

The closed loop singular value plots are shown in figure 4.6, which clearly shows 
that all the plots are well below unity. This indicates successful Tico design. 

Figures 4.7, 4.8 and 4.9 show the time responses (for the first two seconds) for a 
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CLOSED-LOOP MUrFull Order Controller 



Figure 4.3: The HIMAT Closed Loop ji with Pull Order Controller 


Closed Loop MU plots-With reduced order controllers 



Figure 4.4: The HIMAT Closed Loop jj. with Reduced Order(/iA(P)) Controller 
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Figure 4.5; The HIMAT Closed Loop fj, with Reduced Order(^C'(P)) Controller 


Closed Loop SinguUif’ Value Plots for full and reduced order controllers 



Figure 4.6: The Closed Loop Singular VaIues:continuous-full order; dashed-^A(P); 
dotted-/LiC(P) 
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Figure 4.7: Step Response with Pull Order Controller for Nominal System 


Closed Loop Time Response with Reduced Older Controller 



Figure 4.8: Step Response with iJ.A{P) Controller for Nominal System 
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B’igvire 4.13: 'FIk' HIMAT Closed Loop fj. with 8th Order Controller 


step input into the* first channel, for the nominal plant. The responses are very similar 
with slight oscillations in l.lie reduced order plots which are quickly damped out. This 
again indicates f.liat t,he reduction approaches are successful. Time responses for the 


same input signal for the plant perturbed by Aq = 


0.1 0 
0 - 0.1 


is shown in figures 


4.10, 4.11 and 4.12. These figures show that the reduced order controller performs as 
well as the full order controller even for the perturbed plants. 

The corresponding results for the optimized 8*^ order controller are presented 
through the figures 4.13 to 4.16. While the closed loop peak n has shot up to just 
above 1.1, tin; singular value plots are still well below unity. The step responses are 
as good as the full order controller for the nominal as well as the perturbed plant. 
Thus the proposed techniques have successfully achieved a reduction in the controller 
order from 30 to 8 preserving robust stability and performance. 




Ctosed Loop Mawmiim Singular Valua with 8th ordar Controller 



Figure 4.14: The Clo.sed Loop Singular Value with 8th Order Controller 


Nominal Ciosad Loop time response with 8th order controller 



Figure 4.15: Step Response with 8th Order Controller for Nominal System 
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Perturbed Closed Loop Tiine Response wim 8lh Order Controller 



secs 

Figure 4.16: Step Response with 8th Order Controller for Perturbed System 

4.3 Reduced Order Autopilot design for flexible 
Launch Vehicle 

4.3.1 Introduction 

This example consists of two parts. In section 4.3.2 a full order controller for the 
experimental data available for the launch vehicle is designed according to specifica- 
tions. In the second part(section 4.3.3) we apply the proposed reduction algorithms 
to get reduced order controllers. The reduced order controllers thus achieved are 
found to meet all the performance specifications as guaranteed by Chapter 3. 

4.3.2 Plant Description 

The data for the launch vehicle is provided in terms of three (4, B, C, D) matrices 
representing the Nominal Plant {Gnom), the upper perturbed plant((?„) and the lower 
perturbed plant(Gd- Th® upper and lower perturbed plants represent the extreme 
plants and all other plant descriptions axe expected to he between these two hmits. 
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Figure 4.17; The Plant and Controller Interconnection with Weights 


The (tI, B, C, D) matrices incorporate the launch vehicle rigid body, slosh and launch- 
ing modes along with actuators and the nozzle modes. The closed loop representation 
of the Plant and Controller interconnection is shown in figure 4.17. The dashed box 
named G represents the actual system consisting of the nominal model Gnom along 
with the uncertainty weight Wg and the unstructured full block multiplicative un- 
certainty A.Q. The input signals dto 2 ] along with Wg^st represent the wing gust 
disturbance while Send = represent the reference signal. V’per/ and ttper/ 

are the performance variables created for converting the control objectives into TCco 
design objectives. The nominal plant Gnom consists of: 


19 states 
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one control input 


(5-angle of deflection 


two measured outputs 


V'-attitude 
•^-attitude rate 


The control objective is that ipmeas should track the command signal 5cmd- In the 
Tioo design perspective this objective is transformed into a transfer function mini- 
mization problem by suitable selection of weights. To convert the original control 
objectives into the Tioo design objective we use a model following approach with the 
ideal model providing the required i'espoiise{tpideai) of launch vehicle. The final 
design objective is thus minimization of the transfer function from [rjcmd dwi to 
the performance signals [ipperf '‘^perfV- 


Exogenous Disturbance, Commands: There are two sources of exogenous 
signals-wind gusts and angle of deflection commands. 


Wind Gust: The Wind Gust data being unavailable, it is modelled as a direct low 
frequency disturbance on the controlled variable ipmeas- The weights Wgust is 


chosen to be 



0.01(0.5s+l) 

is+l) 

0 


0 

0.01(0.5s+l) 

is+l) 


This implies that the disturbance is expected to be within a bandwidth of 2 


rad/s and the maximum amplitude expected is 1% of the command signal. 
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Deflection Command: The guidance system produces Deflection angle commands 
which the autopilot must make the launch vehicle to follow. Typical deflection 
commands are modelled as ~ ^^crnd'Homd where ?7cmd 1^ assumed to be amy 
arbitrary signal with |l?7cmd|loo ^ 1- Here the Wcmd is chosen as: 


HcTTid — 


0.5(s + 2) 

(5 + 1 ) 


The particular choice roughly implies that the command input is expected 
within a bandwidth of 2 rad/s. This being the bandwidth of the autopilot 
itself, this choice is justified. 


Errors: Two separate variables must be kept small in the presence of the above 
disturbances. In this context these variables will be considered errors. 


Control Input: This variable must be under 5 radians according to specifications. 
So the coiitrol weight is chosen as 


^control — 


1 

5 


The Attitude Error: The ideal attitude response to the angle of deflection com- 
mand is chosen to be 


“^ideal 


s2 -I- 2 X 0.7 X 2s + 4 


^cmd 


Thus our ideal plant has a natural frequency of 2 rad/s and a damping ratio of 
0.7. The tracking error is defined as V'per-/ = i’ideal ~ i’meas- The tracking error 
signal is weighted by frequency depended weight Wp which is chosen as 


Wp = 


0.5(s + if 

(s + 0.1)2 
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This weight emphasizes the need for low tracking error at low frequency. This 
restriction is relaxed slowly for higher frequencies. The actual gain and the 
exact form of the weight is arrived at iteratively and is somewhat arbitrary. 

Uncertainty Model: The given data severely restricts any chance of non- 
conservative modelling of the uncertainty. As haa already been discussed, the data 
is available as a nominal plant with upper and lower perturbed limits. Only way to 
model this sort of uncertainty is as unstructured full block multiplicative type. We 
describe the set of plants by 

Q ’■= {G'nom(l + Ac?Wde/) : stable, IIAgIIoo < 1} 


where Waci is chosen as 


Wdel 


5 ( 5 - 1 - 1 )^ 


( 5 + 10)2 

and Aq is full block type. The Wdei chosen allows for small amount of uncertainty in 
the plant at low frequencies and the amount of allowed uncertainty increases smoothly 
at high frequencies. The particular choice of Wdei ensures that any plant G €0 also 
belongs within the limits that has been set by the upper and lower perturbed plants. 

Now that all the different weights have been specified we present a schematic 
diagram in figure 4.18 of the augmented plant in the generic feedback arrangement. 

The augmented Gaug consists of 28 states with 6 outputs and 5 inputs. The 
controller K is designed using D-K iteration with the “dkit” routine of MATLAB. 
The best Kf found is of 38 order with three inputs and one output. To check robust 
performance we define the augmented uncertainty block according to theorem 1.8.3. 


Aa[/g = 


Ag 0 

0 A. 


: Ag€C1^\ApGC 


'^3x2 
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Figure 4.18: The Augmented Plant with Controller in the generic feedback arrange- 
ment 

The upper and the lower bounds of pL{Ti{GAUGi^f)i^AUG) are plotted in figure 
4.19. The supremum over frequency of the upper bound is found to be 1.082, which 
is acceptable for our purposes. This 38*^ order controller is found satisfactory in all 

respects as will be verified through various time simulations in the next section. But 
such high order controllers are hardly implementable, even on a computer. Thus it is 
a suitable candidate for testing the controller reduction algorithms of Chapter 3. 

4.3.3 Reduced Order Controller design 

It should be noted that both our reduction approaches fiA{P) and nC{P) assumes 
the full order controller satisfying performance specifications to be known. We use 
Kf of previous section for this purpose. The results of the controller reduction are 
presented in table 4.3. 

Figures 4.19, 4.20 and 4.21 show the /x plots for the full order and the reduced 
order controllers. All are found to be in the prescribed range, as had been guaranteed 
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Reduction Approach 

Controller Order 


Pull Order 

38 

1.082 

M(-P) 

16 

1.2072 

/iCCP) 

17 

1.1713 


Table 4.3: Results for Controller Reduction 



Figure 4.19: Closed Loop /r with Full Order Controller 


a priori by proposition 3.2.2 and 3.3.1. The maximum singular value plots for the 
closed loop system with the different controllers are shown in figure 4.22. AH three 
singular value plots lie below unity and thus indicate successful Ttoo design. 

The various controllers in closed loop are now tested in time domain. Step re- 
sponses to the tracking command are shown for the nominal and perturbed closed 
loop with the various full and reduced order controllers. The wind gust disturbance 
is modelled as a step input after 32 seconds. The full order step response is satis- 
factory. It not only gives a well-damped and tolerably quick response with a rise 
time of approximately 2 seconds and setthng time of approximately 4 seconds, it also 
successfully rejects the wind gust disturbance. The peak overshoot is well within the 
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Fr«q (rad/s) 

Figure 4.20; Closed Loop fj, with Reduced Order (/iA(P)) Controller 


Closed Loop mu plots with reduced order oonlroUer 



lO"* 10'* 10‘* lO” 10* 10’ 10* 10* 10* 

. Freq (red/s) 


Figure 4.21: Closed Loop /x with Reduced Order(/xC(P)) Controller 
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Figure 4.22; 
dotted- IJ,C{P) 


Oosed Loop Singular Value Plots 



Closed Loop Singular Valuesicontinuous-full order; dashed- fj.A{P); 


Step Response with Futt Order Controller 



Figure 4.23: Step Response with Full Order Controller for Nominal Model 
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St«p Response with Full Order Controilar 



Figure 4.26: Step Response with Pull Order Controller for Upper Perturbed Model 


Step Response with Reduced Order Controller 



Figure 4.27: Step Response with iJ,A{P) Controller for Upper Perturbed Model 
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Slap Response with Reduced Order Controller 



Figure 4.28: Step Response with ij.C{P) Controller for Upper Perturbed Model 


Step Response with Full Order Controller 



Figure 4.29: Step Response with Pull Order Controller for Lower Perturbed Model 
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Stop Rasponso with Raducad Order Conlrotler 



Figure 4.30: Step Response with fjiA{P) Controller for Lower Perturbed Model 


Slap Rasponsa with Reduced Order Controilar 
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Figure 4.31; Step Response with /iC(P) Controller for Lower Perturbed Model 
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specified limit of 20%. The step responses with the reduced order controllers verify 
that the reduced order controllers perform as good as the full order controller. This 
has already been guaranteed by the corresponding fj, plots and is borne out by the 
time responses. Thus the proposed controller reduction schemes successfully achieve 
considerable reduction in controller order (from 38 to 16 or 17) preserving robust 
performance and stability. 

4.4 Conclusions 

In this chapter the controller reduction techniques proposed in Chapter 3 has been 
tested on two practical examples. The full order design for the experimental HIMAT 
aircraft is taken from [8] and the proposed controller reduction algorithms are tested 
on this available design. The thirtieth order controller is successfully reduced to a 
fourteenth order controller without any optimization and this reduced order controller 
is found to have comparable performance with the full order controller. In the second 
part a complete /x-synthesis design is performed for the experimental data available 
from a flexible launch vehicle. The resulting 38th order controller is found to provide 
the desired characteristics. On application of the reduction algorithms a 15th order 
controller is found that performs comparably with the full order controller. 



Chapter 5 
Conclusions 


In this final chapter we summarize by fisting the contributions of this dissertation 
and recommending the directions for further research in this field. 

5.1 Contributions of this dissertation 

1. In this thesis a new proof has been given for the Kavranoglu’s additive reduction 
technique using the extension of the mixed n theorem for uncertainties with 
frequency dependant bounds of Tits and Balakrishnan [25]. This proof provides 
a rigorous basis for additive controller reduction in the fj, framework. 

2. For the same technique the weights for two sided reduction, which are sufficient 
to guarantee the closed loop ji to remain less than unity, have been derived. It 
is also shown that these weights satisfy the desirable condition that their size 
is minimized in some sense. 

3. One of the main drawbacks of the above method is that it assumes the controller 
to be approximated is stable. Robust controller design techniques like Tioo loop 
shaping or ^.-synthesis do not guarantee stability of the controllers. Thus a 
new coprime factor reduction scheme has been proposed, which do not require 
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stability of the full order contr^^ii^, Knt . , 

ier but guarantees closed loop stability and 

performance with structured urr<:ig ^ , tv, • r 

ertamty. ihe copnme factor perturbations to 

the controller have been found to kn™ a v^l i i 

have a block diagonal structure thus ir 

the reduction algorithm. 


! improving 


4. The above algorithn^ have been ^ 

aircraft. The proposed algorithm^, have been fonnd to work qnite well producing 

nior6 tliStii 50% ^^duction in thf^ Ors j. n 

Controller order without optimization. Detailed 

simuls^tiioii results bsen i *1 

en presented with the full and the reduced order 

controllers for botli tbe nomin/^,1 j 

tae nominal perturbed systems. 

5. A. full jjj synthesis desisrn bas bpfvv. r i i 

ign nas performed on the experimental data avail- 
able from a flexible launch vehiclK, 'T'u j j. *, p , 

vemci^ The details of the design with the selected 

performance weights and the sim^ i p- i.l 

simulation results are presented. The reduction 

e,re also tested on this v. 4 -* i , 

tnis practical example and found to produce con- 

srderable reduction in the controh^r order while preserving performance. 

5.2 Directions of Furtl^er Research 


In the proposed methods, while the freo, j j r n, i 

niie trie rrec^ugney dependance of the structured singular 

M has been utilized fully the spati^ structure of /x have not been utilized. The 

of the frequency dependant se-^; bounding the uncertainties used in Lemif-a 

very restrictive. The extensiop Qf ^be ^ theorem for frequency dependent 

nds on the uncertainty where the bq^^ding set is more general, may lead to less 

conservative techniques of model reduction 
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Appendix 

HIMAT Data 



-0.0226 

-36.6000 

-18.9000 

- 

-32.1000 


0 

0 

A = 

0 

-1.9000 

0.9830 

0 

B = 

-0.4140 

0 


0.0123 

-11.7000 

-2.6300 

0 


-77.8000 

22.4000 


0 

0 

1.0000 

0 


0 

0 



0 

57.3000 

0 

0 


0 

0 





D = 




0 

0 

0 

57.3000 


0 

0 


The design weights are Ws - fg® , Wp = and W„ = 

Additive Reduction Weights 

_ 3.8918s3 + 4.8483e + 0035^ + 1.5088e + 006s + 2.7701e + 003 
“ s3 + 4.4616e + 003s2 + 4.9755e + 006s + 1.8748e + 003 

Multiplicative Reduction Weights 

0 -72378^ + 458.3224s2 + 7.2567e + 004s + 169.8140 
luvis) - ^3 + 4.i844e + 003s2 + 4.3773e + 006s + 3.9744e + 004 
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Launch Vehicle Data 
Additive Reduction Weights 

9.6768s 3 ^ 2.5500e + 0043^ + 9.4549e + 005s + 8.5589e + 006 
" s3 + 1.6302e + 003s2 + 7.4495e + 005s + 1.0127e + 008 

Multiplicative Reduction Weights 

0.0866s 3 + 5.0172e + OOOs^ + 1.8747e + 005s + 1.0746e + 006 
7t/v(s) - ^3 3.6585e + 003s2 + 8.2348e + 005s + 3.8928e + 007 




